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Abstract 

The three-phase version of the hybrid chiral bag model, containing the phase of asymp- 
totic freedom, the hadronization phase as well as the intermediate phase of constituent 
quarks, is proposed. For this model the self-consistent solution, which takes into ac- 
count the fermion vacuum polarization effects, is found in (1+1) D. Within this solution 
the total energy of the bag, including the one-loop contribution from the Dirac's sea, is 
studied as a function of bag geometry under condition of nonvanishing boson conden- 
sate density in the interior region. The existence and uniqueness of the ground state 
bag configuration, which minimizes the total energy and contains all the three phases, 
are shown. 
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1 Introduction 



The idea to describe a hadron as a bounded region of space ("the bag ") filled with quark 
and gluon fields appears nowadays to be one of the most natural ways of constructing the 
effective quantum field theory with total confinement of color objects [1-12]. The boundary 
conditions are chosen as to confine the colored particles to the bag without breaking the rel- 
ativistic covariance [1-3]. The very possibility of existence of such an object is related to the 
change of the vacuum structure inside the hadron. It is supposed that the non-perturbative 
vacuum inside strongly interacting particle is completely destroyed, what produces an excess 
of energy proportional to the bag volume [1-5]. The stability of this system is provided by the 
valence quarks contribution to the total energy proportional to the inverse bag radius. Even 
the early MIT-bag model yielded consistent results for the mass spectrum and other static 
characteristics of hadrons with relatively small set of free parameters [6-8]. Further develop- 
ment of this model has led to taking into account the effects of spontaneous chiral symmetry 
breaking — one of the most important features of strong interactions at low energies, and 
to incorporating the meson fields into the theory which play the role of Goldstone bosons 
[9-11]. (Depending on the number of light-quark flavors, it is either the triplet of 7r-mesons 
for SU(2) group or the octet of ix- and .fT-mesons for SU(3)). As a result the most consistent 
approach to the description of hadronic structure based on bag models has been elaborated 
within the framework of so-called Hybrid Chiral Models (HCM) [10-12]. The HCM treat the 
nucleon as the little bag with quarks and gluons confined inside, surrounded by a big cloud 
of virtual mesons. The latter can be described either by the "hedgehog" 7r-meson solution 
of the Skyrme model, or by some more complicated models including vector mesons [13-15]. 

On the whole, the models of such a type predict correct scales for various static char- 
acteristics of non-strange baryons. At the same time, the estimation of certain values may 
yield essential quantitative discrepancies, even the sign may be incorrect [12,16]. Moreover, 
the best HCM results have been obtained for such values of model free parameters (which 
are the vacuum pressure constant B, the coupling constant a s , the current quark masses, 
and the constant of Casimir energy Z), that differ from their values obtained by means of 
other methods [4,5]. These circumstances indicate that the current formulation of HCM 
suffers from a number of shortcomings. 

Actually, the most subtle point in the HCM is the Cheshire Cat Principle (CCP) which 
is basic for parting space into regions with different phases inside of them [17]. The essence 
of CCP is the hypothesis [18] that the fermion theory inside of bag and the boson theory 
outside are actually equivalent and can be transformed into each other via the bosonization 
procedure. Therefore none of physical properties of such bag depends actually on the choice 
of boundary surface, while the boundary conditions are determined from the bosonization 
equations [19]. However, the bosonization, as well as CCP, can be rigorously proved only 
in (1+1)D, while in the real (3+l)D world the solution for the bosonization problem is 
still absent. As a result, in the (3+l)-dimensional HCM based on CCP there exists a rather 
small set of observables (e. g., the topological charge), which do not really depend on the bag 
radius [20] . Moreover, the phenomenology of strong interactions predicts unambiguously the 
existence of the characteristic confinement scale about 0.5 fm, and so in the realistic (3+l)D 
models the CCP should be strongly violated, regardless on the proof of bosonization. 

Therefore it seems to be natural to modify HCM as to avoid any relation to CCP and 
specific features associated with, such as an infinitely thin boundary surface between dif- 
ferent phases and corresponding surface action. Such a modification appears to be possible 
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provided the correlation between the different phases of bag, which are now not assumed to 
be equivalent in the sense of bosonization, arises due to some interaction of the real status 
taking place in a finite region of space. The latter can be naturally treated as the region of 
the intermediate third phase of the bag [21]. The introduction of such an additional phase is 
useful also because it allows to introduce the chirally invariant mechanism of the quark mass 
generation, hence it can be considered to a certain extent as the phase of constituent quarks, 
while the original 2-phase model does not take into account their existence at all. Mean- 
while, the concept of constituent quarks has been shown to be very efficient in the hadron 
spectroscopy. From this point of view the most attractive physically situation is the one, 
where the initially free, almost massless quarks (corresponding to large momentum transfer) 
transmute firstly into "dressed" due to interaction massive constituent quarks carrying the 
same quantum numbers of color, flavor, and spin, and only afterwards there emerges the 
purely mesonic colorless phase. 

The first order approximation to such version of the bag is given by the 3-phase hybrid 
model with additional constituent quark phase instead of the boundary with vanishing radial 
extent [21]. This model allows to take into account three phases: the phase of asymptotic 
freedom with free massless quarks, the phase of constituent quarks which acquire an effective 
mass due to the chirally invariant interaction with the boson field in the intermediate region 
of finite size, and the hadronization phase where the creation of free quarks is suppressed by 
infinitely large mass, while the non-linear dynamics of boson field causes the appearance of 
the c-number boson condensate in the form of a classical soliton solution, which ultimately 
accounts for the quantum numbers of the whole bag. 

In the present paper the toy model of such kind in (1+1)D is considered, where a single- 
flavor fermion field is coupled in a chirally invariant way to the real scalar field possessing a 
non-linear soliton solution in the exterior region. The self-consistent solution to the model 
equations, which takes into account the effects of fermion-vacuum polarization, is found. 
Within this solution the renormalized total energy of the bag is studied as a function of the 
bag geometry under condition of nonvanishing boson field in the interior region. It is shown 
that by suitable choice of the model parameters the configuration with minimal total bag 
energy and containing all three phases, exists and is unique. 

2 Lagrangian and Equations of Motion 

The division of space into phases is performed by means of the system of subsidiary fields 
9(x) according to the method of ref. [3, 22]. The essence of this method may be explained 
as following. Let us consider the Lagrangian of the form: 

C = \ [d^f - 9V{<t>) + \ (d,9) 2 - g 2 W(9), (2.1) 

where the coupling constant go of self-interaction of field 6 is assumed to be large enough 
to neglect the matter fields in the dynamics of 9 to the leading order, and thereafter to 
use 9 as background fields for the dynamics of 0's [21,22]. One can obviously construct the 
Lagrangian containing as many fields 9{x) as needed with appropriate self-interaction, which 
will determine (almost) rectangular division of space into regions corresponding to different 
phases, while the Lorentz-covariance will be broken only spontaneously, namely on the level 
of solutions of equations of motion. The latter circumstance allows one to use the framework 
of co variant group variables [23] in order to restore the covariance. Assuming further that 
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subsidiary fields 9{x) have already formed the required bag configuration, let us start with 
the following Lagrangian: 

C = $idrf> + \ {d^f - 9{ Xl < \x\ < x 2 ) [ij, e i9 ^ij] - 

~0(\x\ > x 2 ) [i>,e l9 ^}_ + V(<p)j . (2.2) 

The commutator of fermion fields in terms, describing the chiral fermion-boson coupling, 
ensures the charge conjugation symmetry of the model. 

So in initial stage we have the theory of two fields, the spinor field ip and the boson 
field (p. These fields are free and massless in the region I : {\x\ < Xi}; in the region 
II : {xi < \x\ < x 2 } the field <p interacts in a chirally invariant way with fermions resulting 
in emergence of the effective fermion mass M; and in the region III : {\x\ > x 2 } the effective 
fermion mass increases up to Mo and the self-interaction of ip switches on, providing the 
appearance of a soliton solution for the boson field in this region. Note that in this model the 
vacuum pressure term appears to be redundant, since due to the existence of the intermediate 
phase the Dirac's sea polarization behaves very specifically and ensures itself the required 
"inward pressure". Moreover, in our model there is no special need in "valence" fermions, 
since it is the boson condensate in the form of the topological soliton, which accounts now 
for all the bag's quantum numbers, represented by the topological "baryon number" in this 
simplest case. The latter point makes it difficult to explain the need for the vacuum pressure 
term in this model too. 

In order to provide the confinement of fermions we assume the mass M to be very large, 
what leads to the dynamical suppressing the fermion field in the exterior region At 
the same time, in the region of asymptotic freedom we have free (decoupled from fermions) 
massless scalar field and so the nonvanishing density of the boson condensate. This situation 
does not contradict the general concept of bag models, and can be treated as one of possible 
versions of the 3-phase model. Another possibility with vanishing scalar field in inner region 
has been studied in detail in ref. [21]. Let us assume also, that the solution of equations of 
motion for the boson field is an odd topological soliton. The case of an even soliton with 
nonvanishing scalar field inside is of no interest, since in this case the trivial solution (scalar 
field is a constant equal to one of minima of self- interaction potential V(ip)) is the only 
energetically preferable. 

Let us consider now the behavior of fields in detail. According to the general approach 
accepted in hybrid models we consider the boson field in the mean-field approximation, i. e. 
it is assumed to be a c-number field. Neglecting temporarily the explicitly Lorentz-covariant 
description, we'll consider the center-of-mass system of the bag, where (p(x) is stationary 
classical field being a background for evolution of fermions. The equations of motion read: 
in the region / 

ibi) = 0, (2.3a) 
<p" = 0, (2.36) 

in the region II 

(id - Me i91hip ) i> = 0, (2.4a) 



M 



4>, l5 e i9 ™4> ), (2.46) 
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and in the region III 

(id - M e i915<p )if; = 0, (2.5a) 

-</ + V'((p) = , (2.56) 

where ( ) in eq. (2.4b) stands for the expectation value in a given bag state. To simplify 
calculations, we put further g — 1, because the dependence on it can be easily restored by 
means of the substitution (p — > tp/g. Then the spectral problem for fermionic wave functions 
ipui with definite energy u reads 

ufa = -ionj}' u + (3e?^ [M6( Xl < \x\ < x 2 ) + M 6(\x\ > x 2 )\ (2.6) 

To complete the formulation of the spectral problem, the following boundary conditions have 
to be imposed 

±i 1 l ^(±x 2 ) + e^^ty^ixa) = , (2.7) 

provided with the condition of continuity for ip(x) on boundaries separating the regions / 
and II. Note that the boundary conditions (2.7) are actually the standard chiral boundary 
conditions for the hybrid models [9-12]. However, they arise now as a direct consequence 
of an infinite mass of fermions in the region rather than from the local surface action, 
what is not completely correct [21]. In the region I the equation (2.6) is the equation for 
free massless fermions 

ufa = -iaipj , (2.8) 
while in the intermediate region // one has 

oo^jj = -iutfu + fiMe^ipn . (2.9) 

Conditions of wavefunction's continuity on the boundary between I and II are given by 

while at points | X | — X2 the wavefunctions satisfy the conditions (2.7). By this the field 
if in eq. (2.9) is not arbitrary but has to be determined self-consistently from eq. (2.4b) 
with appropriate continuity conditions imposed on the field and it's derivatives at points 
\x\ = x 1>2 . 



3 Self-consistent solution of bag equations 

The essential feature of this bag configuration is the fact, that the self-consistent equations 
(2.4) in the closed intermediate region II of finite size d = x 2 — X\ possess simple and phys- 
ically meaningful solution, which would be unacceptable if these equations were considered 
in the infinite space. In order to obtain this solution, we perform in the region II the chiral 
Skyrme rotation 

ip = exp(-i7 5 vV2)x, (3-1) 
hence the equation (2.9) and the boundary condition (2.7) transform into 

(u - \ <pf) Xu = -iaxL + PM X „, (3.2) 
±i 1 l Xu {±x 2 ) + Xu{±x 2 ) = (3.3) 
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correspondingly. 

It follows from (3.2) and (3.3) that if we assume the linear behavior of the scalar field in 
the region //, namely 

if' = const = 2A, (3.4) 
then the equation (3.2) becomes the equation for free massive fermions: 

vx = -iocx' + pM X (3.5) 

with eigenvalues v = us — \. Therefore, as we expected, the massless (in the region I) fermions 
acquire the mass M in the region 77 due to the coupling to field ip, and so the intermediate 
phase emerges which describes massive quasifree "constituent quarks" . Continuity conditions 
for the scalar field, which should be an odd function, lead us to the unique solution 

ip(x) = 2Xx . (3.6) 

The equation (3.5) obviously possesses the sign symmetry v — > —u, which corresponds to 
the unitary transformation of fermionic wave function 

X -»• X = hiX • (3.7) 
The important point here is, that the axial currents 

h = #75e n5 ^ = ix + liX (3-8) 
coincide for these sign-symmetric states 

h = iX + liX = ix + liX = h ■ (3-9) 

In general, however, one cannot derive the sign symmetry for the fermion spectrum in our 
case from the sign symmetry v <-> —v in eq. (3.5), because the latter is true only in the region 
//, while the spectrum has to be determined from the Dirac equation on the unification of the 
regions / + II. The straightforward solution of equations (2.8-9) with account of boundary 
conditions (2.7) and constraint (2.10) gives the following equation for spectrum: 

y _ c -2ikd M-i(v+k) j _ c 2ikd M-i{v~k) 
I a- r> • \ ' M—itu—k) M—i(v+k) /0 , A \ 

exp(4^ 1 - 2 m ) = M+l(u . k) ^ kd M +t( u + k) > ( 3 - 10 ) 

x c M+i(u+k) x c M+i{v-k) 

where v 2 = k 2 + M 2 and </?i = ip(xi). Analysing the equation (3.10) one easily finds that 
the fermionic spectrum reveals the symmetry v <-> —v, if 

A\x\ — 2(fii = tts , (3.11) 

where s is integer, since for such values of the derivative of the field (p(x) the 1. h. s. 
of equation (3.10) reduces to (— l) s exp(4iz/xi). However, in our case ipi = 2\x±, hence 
the equation (3.11) leads to the single possibility s = 0, while the parameter A remains 
arbitrary. This result differs crucially from the case considered in [21], where tpi — by 
virtue of vanishing of the boson condensate in the region J, and so the set of solutions 
with different s ^ emerges. In the latter case, however, the equation (3.11) produces the 
non-trivial relation between A and x±. 



6 



According to eq. (2.4b), in the region II (p"(x) is determined by the v.e.v. of the C-odd 
axial current 

.h = \ [Vv 75 e l75 V]_ = | [x + ,*7ix]_ , (3.12) 
with x being now the secondary-quantized Dirac field in chiral representation (3.1) 

X (x,t)=Y / bnXn(x)e-^ t , (3.13) 

n 

where Xn( x ) are the normalized solutions of the corresponding Dirac equation, and b n , 6+ are 
fermionic creation-annihilation operators, which obey the canonical commutation relations 

{b n ,b+,} + = <W , {b n ,b n >} + = . (3.14) 

The average over the given bag's state includes, by definition, the average over the filled sea 
of negative energy states u n < + possible filled valence fermion states with u n > , which 
are dropped for the moment because their status is discussed specially below. Finally, 

(J 5 > = (J 5 )sea = [\ E -\ E ) Xj<7lXn " ( 3 ' 15 ) 

Let us emphasize, that in (3.15) the division of fermions into sea and valence ones is made in 
correspondence with the sign of their eigen-frequencies uj n , which differ from sign- symmetric 
v n by the shift in A 

co n = u n + X , (3.16) 

and so do not possess the symmetry uo <-> —uo. However, if we suppose additionally that 
u n and A are such that for all n the signs of v n and uj n coincide, i. e. after shifting by A 
none of z/ n 's changes its sign, then the condition w n ^0 in eq. (3.15) will be equivalent to the 
condition u„ ^ . Hence 



< 



(J 5 )sea= U E ~\ E tfi7lXn = (3.17) 

by virtue of relation (3.9). In turn, it means that the eq. (2.4b) in the region II reduces 
to ip" = 0, what is in excellent agreement with the assumption that tp'(x) = const in the 
region //. In other words, we obtain the solution of coupled equations (2.4) in the region // 
in the form of the linear function (3.6) for the scalar field, and the equation (3.16) for the 
fermion energy spectrum, where v n is defined from eq. (3.10) after replacing the 1. h. s. to 
exp(4iz/xi). 

There are the following keypoints, that make this solution meaningful. The first is the 
finiteness of the intermediate region size d, because for an infinite region // the solution (3.6) 
is obviously unacceptable. In our case, however, the size of the intermediate region is always 
finite by construction, and the boson field tp(x) acquires the solitonic behavior in the region 
/// due to the self-interaction V(ip). Here the following circumstance manifests again: in 
(1+1)D the chiral coupling ipe 115 ^^ itself cannot cause the solitonic behavior of the scalar 
field by virtue of the effects of fermion-vacuum polarization only, i.e. without additional 
self- interaction of bosons [25]. 

Another important circumstance is the calculation of the Dirac sea average of axial 
current J5 by means of the symmetry v <-> —v. The point is that for such kind of averages 
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another consistent definition is possible, namely, via so called 77-invariant, which is often used 
for the study of fermion vacuum polarization effects [26]: 



lim 



-|w„|?7 + 



(3.18) 



Since the spectrum of u n is not symmetric, the expression (3.18) doesn't vanish, in contrast 
to the average (3.17). The expression (3.18) with connection between u n and u n taken into 
account, transforms into 



lim 

77^0 



sinh A77 ^ 



e ^XnhlXn 



(3.19) 



and for r\ — > the sum in (3.19) diverges as I/77, so such (J 5 ) appears to be proportional to 
A and doesn't vanish unless A 7^ 0. The same makes the r.h.s.of eq. (2.4b), and so the linear 
functions will not be the solution of this equation in the region II. 

Moreover, this is true for other fermionic averages too, in particular for the v. e. v. of 
the fermion charge. In our approach, by virtue of discreteness of the fermionic spectrum, 
the vacuum value of the charge is zero, in analogy to the case of the axial current. Namely, 
if the C-odd expression for the charge 



(3.20) 



is used, then under the same conditions of sign-correspondence between oj n and v n for each 
n, one obtains for the average of Q the following 



(Q) = {Q)t 



SE-SE=SE-SE=o- 



(3.21) 



UJ n <0 OJ n >0 Vn<0 Vn>0 

However, if (Q) sea is defined via 77-invariant, then one obtains 



{Q)sea = lim 

•q— >U 



5 £• 

Wn<0 



<^n>0 



7^0 



(3.22) 



for A 7^ by virtue of absence of the symmetry u <-> —uo. Nethertheless, there are serious 
reasons to consider the relations (3.17) and (3.21) as the most adequate way of calculation 
the sea average in our problem. First, 77-invariant is actually the measure of the Hamil- 
tonian asymmetry, rather than the fermion charge of the ground state. Second, the other 
regularization schemes besides of thermal regularization are possible for divergent sums such 
as (J5) and (Q). Also, there is no warranty that the regularized expressions (3.18) and 
(3.22) would produce the correct results for 77 — > 0. It is easy to give an example, when 
the dependence on a parameter in the sum or integral is not continuous, in particular, 
F(k) = f£°dx sinkx/x = ix/2 for any k > 0, but F(0) = 0, i.e. in this case the limit of the 
integral for k — > and it's exact value at k — are different. 

Let us consider further an adiabatic process of changing the gradient of the boson field A 
in the region 77, starting from A = 0. In the initial moment uo n = is n , so any C-odd fermionic 
sea average obviously disappears: (. . . ) sea = for any reasonable regularization scheme. 
Now let us take into account that in our case the fermionic spectrum is discrete and depends 
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on A continuously, so for A small enough all the u; n 's keep their signs. Then from general 
grounds (. . . ) sea should vanish still, although the spectrum u n is already non-symmetric 
and so the 77-invariant is nonzero. Therefore it is reasonable to assume that in this case 
the thermal regularization can yield physically inconsistent results for C-odd quantities such 
as the axial current and the fermion charge of the vacuum. It is worth-while to note the 
essential difference of this situation from the case of infinite space, when the spectrum is 
continuous. In the latter case, the density of states changes for any small change of A, and 
so all the averages change too, what gives rise to the effect of induced fermion numbers [26]. 

Finally, it is possible to check the self-consistency of solution (3.6) numerically. In this 
case the problem is solved on a lattice, and so the number of degrees of freedom, as well as 
the number of fermion levels, is finite, so one doesn't need any regularization at all. The 
results of these calculations show that the linear function (3.6) and the vanishing averages 
are the only self-consistent solution to this problem. 

The interpretation of this result is as following. In our case the boson field has no 
discontinuities at all (it doesn't vanish in the region /, too), and is topologically equivalent 
to the odd soliton which would take place in absence of fermions due to the self-interaction 
V(ip) only. Therefore, the topological number does not depend on the existence and sizes 
of the spatial regions with fermions (/ and //). On the other hand, the baryon charge of 
the hybrid bag is, by definition, the sum of the topological charge of the boson soliton and 
the fermion charge of the bag interior. The latter one is zero in our case, hence the baryon 
charge is determined completely by the topological charge of the boson field, and so does 
not depend on the sizes of the regions I and //, in correspondence with the general ideology 
of hybrid models. Therefore, in our approach the "hadron" is a composite particle, which 
consists of the boson soliton, to which the fermion bag is coupled through the interaction in 
the intermediate region II. Note also, that although quantum numbers of such a composite 
particle are completely determined by the soliton, it doesn't mean that the filled fermion 
levels with positive energy should not exist at all. This could take place for small enough 
values of the parameter A only. If A increases, the negative levels u) n = — \u n \ + A will 
unavoidably move into the positive part of the spectrum. The change of sign of each such 
level will decrease (Q) sea by one unit of charge, but if we fill the emerging positive level with 
the valence fermion, then the sum Q va i + Q sea remains unchanged. Analogously, the total 
axial current will be equal to J va \ + J sea and will not change too, what ensures the vanishing 
r.h.s. of (2.4b) and so preserves the status of linear function (3.6) as the self-consistent 
solution of the field equations. In other words, by definition the ground state of the bag is a 
state in which all the levels with v n < are filled (the inequality is strong, because there are 
no levels with v = in the considered case). The existence or absence of valence fermions in 
such construction of the ground state of the bag depends actually on the relation between A 
and \v n \min, and so appears to be a dynamical quantity like the other parameters of the bag 
(the size and mass), which are determined from the total energy minimization procedure. 

4 The Total Energy of the Bag 

On the unification I + II the boson soliton takes a form of the linear function (3.6). This 
function (after rescaling (p — > ip/g) is sewn together with the soliton solution of eq. (2.5b) in 
the exterior region by means of the continuity conditions for the field and its derivative. To 
preserve the generality of consideration, in the region /// we'll use the asymptotic expansion 
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of the soliton solution of eq. (2.5b) for large \x\, namely 

Vsol (x) = ^(l-Ae- mx ) , x>x 2 , (4.1) 

with m being the meson mass in the exterior region of the bag, while for x < —x 2 ip so i(x) 
is determined by oddness. The factor 7r/g means that we deal actually with a phase soliton 
with topological charge being multiple of 2ir/g, since it is the period of the initial chiral 
interaction ip exp(i / y 5 g(p)ilj. The constant A is determined from the continuity condition for 
boson field at points x = ±x 2 . Note that in the exterior region there is no chiral invariance 
due to the phenomenology of strong interactions on the one hand, and on the other due to 
the specific features of (1+1)D scalar models, which make the presence of the meson mass 
the necessary condition for the required soliton profile to be formed. 
The continuity conditions at x — ±x 2 give 

2Xx 2 = 7T (l - Ae~ mX2 ) (4.2a) 

2A = nmAe' rnX2 . (4.26) 
Then one finds the relation between the parameter A and the bag's size x 2 : 

2A = tt — , 4.3 

mx 2 + 1 

whence the total energy of boson soliton can be represented as 

E v = (4.4) 
g 2 mx 2 + 1 

The total energy of the bag consists of E v and of the fermionic contribution 

Ebag = E v + E^,. (4.5) 

It is obvious from (4.4), that the energy of the boson field is a smooth function, which 
decreases for x 2 — > oo, producing no vacuum pressure in spite of the fact, that the gradient 
of ip in the domain I + II yields the constant positive contribution to the energy density 
\ (p' 2 = 2A 2 , which could be identified with the vacuum pressure B in the standard HCM. 
Actually, it is an artifact of one spatial dimension in our problem: when the bag's size 
increases, the gradient of ip in I + II will always decrease in any number of space dimensions 
by virtue of (4.3), while the volume of domain I + II in l(space)D increases only linearly 
and so cannot compensate the decreasing of A, what takes place in 2- and 3-(space)D. So 
in (1+1)D the non-trivial dependence of the total energy Eb ag of the bag on the model 
parameters could originate only from the fermion contribution which is given by the 
sum of the Dirac's sea of filled negative energy states and positive energy valence fermions 

E^ = E val + E sea . (4.6) 

For the ground state of the bag described above, the sum (4.6) can be reduced to the single 
universal expression by taking into account, that the charge conjugation symmetry dictates 
the following definition of the Dirac's sea energy [25,27] 

E sea = \ J2 W « ~ I ( 4 - 7 ) 

U>n<0 UJn>0 
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If the transition from uo n to v n is sigh-preserving for all n, and so there are no valence fermions 
in the ground state of the bag (in order to provide the vanishing v.e.v. for the charge and 
axial current), one finds from (4.7) 

= E sea = \ (~Vn + A) - \ £ K + A) = - E v n ■ (4.8) 

U n >0 V„>0 U n >0 

If the parameter A appears to be large enough, so that the initially negative level u n = 
— \v n \ + A changes its sign and turns into the filled valence state, it is convenient to calculate 
in two steps. First, we consider the contribution from all states with \v n \ > A to E sea , 
which in analogy to (4.8) reads 

Kea = (4-9) 

^m>A 

To this expression the energy of the valence fermions E vai = —\v n \ + A and the contribution 
of the positive levels with u n = ±|z/„| + A should be added, what yields 

= -\u n \ + A - \ \{-\v n \ + A) + (K| + A)] + E' sea = - v n , (4.10) 

iy n >0 

i. e. the same eq. (4.8) as we have got for the energy of fermions without filled valence 
states. 

It is convenient to introduce a set of new parameters, in terms of which the total energy 
of the bag will be expressed in the most appropriate form. First of all, we introduce the 
dimensionless quantities 

a = 2Mx!, B = 2Md, p = 2Mx 2 , (4.11) 

and consider the eq. (3.10) in their terms. This equation has two branches of roots. The 
first one corresponds to real k and in terms of parameters a and j3 is determined from 

/ /- -\ x x cos 3x + sin Bx , , „ _ 

tan I a VI + x 2 ) = - — — , 4.12 

v 1 y/x 2 + 1 1-cosBx + x sm fix 

where the unknown quantity is the dimensionless x defined through k = Mx, v = My/1 + x 2 . 
The real roots x n belong to the half-axis < x n < oo, since the fermionic wave functions are 
actually the standing waves in a finite spatial box with degeneracy in the sign of k, while 
the corresponding frequencies v n lye in the interval M < v n < oo. 

The second branch corresponds to imaginary k = iMx, v = My/1 - x 2 , < x < 1, and 
can be derived from (4.12) by means of the analytical continuation 



/ / -\ x x cosh Bx + sinh Bx 

tan (ay/l-x 2 ) = - — - - H . 4.13 

V 1 — x 2 cosh Bx + x sinh px — 1 



For this branch < v n < M. 

Therefore, v n and E^ appear to be the functions of two independent dimensionless pa- 
rameters a and (3, the sum of which defines the dimensionless size of the confinement domain 
P 

a + B = p. (4.14) 

Proceeding further, it is convenient to extract the mass of a "constituent quark" M from 
the sea energy and fermion frequencies as a dimensional factor: 

v n = Me n , e n = ^fl + x 2 n , (4.15) 
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hence = —MJ2 n £ n- Introducing further the dimensionless ratio of the two model mass 
parameters 

p = m/2M , (4.16) 
and the dimensionless total energy £b ag = Eb ag /M, for the latter one finds 

£ bag = £4a,f3) + ^^— . (4.17) 
g pp + i 

Let us note, that in this case the independent model parameters are a, (3, which define the 
sizes of the inner bag's regions, while p is determined from (4.14). In this point the present 
version of 3-phase model differs from the one considered in [21], where a, (3 are defined 
unambiguously by p and p. So the total bag's energy depends ultimately on the following 
dimensionless parameters: p, g and a, (3, where the parameter p is fixed by the ratio of the 
masses m and M, while the optimal values of a, (3 for the ground state should be found 
from the condition of minimum of the total energy for given p. 

The dimensionless fermion sea energy which obviously diverges in the upper limit, 
should be regularized, so a consistent renormalization procedure is required. First, let us 
consider the asymptotics of roots of eq. (4.12) in the UV-domain, when x n >> 1. It is 
convenient to rewrite the eq. (4.12) as 

sin aVl + x 2 = 

= \ (Vl + x 2 + x) sin (aVl + x 2 + f3x + 7) + \ (Vl + x 2 — x) sin (aVi + x 2 — (3x — 7) , 

(4.18) 

where 7 = arctanx . One can show that 

Xn(aJ3) = -ft ± ™ + MTIJE l(j/2 + ™)a/p] + W2 + 
p 7T/2 + 7rn 

In the expression (4.19) the first term leads to the quadratic and linear divergences in J2n £ n, 
and the second one contains the logarithmic divergence, while the term with the sine does not 
yield any divergence at all. Therefore, the UV-regularization consists in the compensation of 
the first term and the divergent part of the second term in the asymptotic expression (4.19). 
As a first step we use that obvious fact, that only the difference between two energies is 
physically meaningful, rather than the energy itself. It could seem, that the most natural way 
is to choose the energy of the free fermions sea in the similar "volume" p as the reference point 
for S^. However, it has been shown in [21], that this subtraction fails in our case because 
this energy appears to be infinitely larger than any configuration with (3 7^ 0, even after any 
possible counterterms are added. Such an infinite energy barrier between the regularized 
£f(a,/3) and the free fermion sea corresponds perfectly well to the intuitive feeling that the 
free fermions can hardly be a reasonable first approximation to the confinement problem. 

As a result, in this case there is no unique prescription for the choice of subtraction point 
in the renormalization of £^(a,(3), what is actually the common feature for the majority 
of bag models [10,12,28]. In the "classical" renormalization procedure, the uncertainty in 
the choice of subtraction point is cancelled by fixing the physical values for corresponding 
number of parameters. For obvious reasons, we won't do that in our "toy" (1+1)D model, but 
instead will consider the most straightforward approach to the compensation of divergences 
in the sum (4.10), which preserves the continuous dependence of the result of subtraction 
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on the model parameters. This approach is based on the subtraction from J2 n £ n the other 
sum with the same summation index n, the common term of which coincides exactly with 
the divergent part of asymptotics (4.19), resulting in the finite quantity 



£</> - 



'tt/2 + nn 1 + 13/2 



7r/2 + Tin / 



(4.20) 



This method requires no counterterms because all the divergencies are already cancelled by 
the subtracted sum. Of course, to some extent the physical meaning of such procedure is 
lost. However, it should be emphasized, that it is only the (1+1)D case, when the theory 
with coupling Cj = Gift (a +1^577)1/1 is (super)renormalizable and any counterterm has explicit 
physical meaning. For higher space dimensions this is already not true, and so the procedure 
of compensation of divergences in the energy based on (4.20) should not be considered as 
having no motivation. For more detailed discussion on the extraction of finite part from the 
divergent Dirac's sea energy see refs. [28-31]. 

Proceeding further, let us turn to the study of the total bag's energy 



£ bag = £*(a, (3) + ^ (4.21) 
g 2 fip+1 



as a function of parameters a, 13. The analysis of the contribution of convergent logarithmic 
part from the sine-term in the asymptotic expression (4.19) to £</, yields the first feature of 
Sbag- Let us transform this contribution to the form 

{S<\ (a,^) = lE(-ir Sin[(WP)(r ! + 1/2)] (4-22) 

V " ' 71+ 1/2 ^ ' 

and then use the wellknown relation 

£(-l)" SiD[ : ( ;t / f )1 - 1 °^/4 + ./4). (4.23) 



n=0 



It is easy to see, that the sums (4.22) and (4.23) have the similar common term, and the 
sum (4.23) diverges as — ln(7r — z) when z — > n. Hence, for na/p — > n , what means either 
(3 — > 0, or a — > 00 for finite /3, the sum (4.22) will show the similar behavior, namely 

(£A (a, (3) -> --\n((3/a) , (3/a -> 0. (4.24) 

V / log 7f 

Therefore, both the regularized fermion energy (4.20) and the total bag's energy possess 
the logarithmic singularity for f3 — ■> 0, and at the same time the logarithmic increase for 
a — > 00 and finite f3. Note that although the possibility of smooth transition to a 2- 
phase configuration, when d — > 0, exists formally for the initial Lagrangean (2.2), actually it 
disappears due to the singularity for (3 — > 0. In other words, in such 3-phase model the radial 
extent of the intermediate phase can be sufficiently small, but not zero, what corresponds to 
the general physical treatment of the structure of many-phase systems. 

£bag will also grow for f3 — > 00 and finite a. Because the increase of /3 yields the increase 
of p, in this case na/p — > 0, hence the logarithmic term (4.22) becomes negligibly small, 
and the main contribution is provided by the next terms in expansion in 1/n. However, it is 
technically more convenient to use the fact that for p — > 00 the fermionic spectrum becomes 
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quasicontinuous with the exception of a small vicinity of Fermi surface with zero energy, 
what allows to transform the sums over x n into integrals over dx. In particular, the analysis 
of distribution of the roots of eq. (4.12) shows, that for this limit J2 n £ n is approximated by 
the following (divergent) integral 



x 



(3 + -—; + a- 



Ve n ^- dx Vl + x 2 x 

7T J 

x 2 sin (a\Jl + x 2 ^j cos (a\J\ + x 2 ^j 



(4.25) 



1 + x 2 x 2 + S i n 2 (a^TTx 2 ) VI + x 2 (x 2 + sin 2 (as/l + x 2 )) 
For the subtracted sum in (4.20) one can easily find 

The integrals (4.25) and (4.26) have the same divergent part 

- [dx {px + 1 fx + /3/2x) , 

7T J 

so their difference yields a converging integral, in agreement with the subtraction procedure. 
The leading term of the intergrand in this difference, taken with the (correct) inverse sign, is 
(3/8jtx 3 . This finally leads to the emergence of the positive, proportional to f3, contribution 
to and correspondingly to Ebag- 

Finally, for a — > and finite f3 there is no singularity in £^ and Ebag at all, what is a 
direct consequence of the analysis performed in ref. [21], where it has been shown that the 
renormalized £^'s for a = and finite a ^ differ always by a finite value. Therefore it 
is the soliton energy and so the current values of parameters p and g, which play the 
main role in this domain. More exactly, for p being not too small, the dependence of on 
g will be such that for sufficiently small g the sharp enough increase of Eb ag takes place for 
a — > and finite /3, providing the existence of a well-defined minimum. On the contrary, 
the sufficient increase of g (or p — > 0) could make to be negligibly small (and almost 
constant) for any bag's size, and then the minimum in the total energy disappears. 

The numerical calculation confirms completely such qualitative predictions for the behav- 
ior of and Sbag- However, due to the absence of singularity in £ bag for a = the question 
about the existence of the minimum of the total bag's energy considered as a function of 
a, (3, and therefore the existence of a stable ground state of the bag itself for given p and 
g can be answered only numerically. Such calculation has been performed for p = 0.25, 
what corresponds approximately to the ratio m 7r /2mg, where the constituent quark mass 
is assumed to be equal to 300 MeV, and for g — 1. The qualitative behavior of £b ag as a 
function of a,f3 is shown on the Fig.l, from which one can easily see, that the total energy 
reveals the unique minimum for non-zero values of a, (3. On the Fig. 2 the equal energy 
curves are shown, what allows to observe this minimum of energy more explicitly. In the 
case under consideration the values of a and (5 in the minimum differ by two orders, i. e. 
the size of the intermediate region appears to be essentially smaller than that of the inner 
one. Therefore to some extent it can be treated as a smeared boundary between the phase 
of asymptotic freedom and the purely colorless phase. On the Fig. 3 the scheme of fermionic 
levels in the vicinity of zero energy is presented for this configuration, which shows that in 
this case the ground state contains one filled valence level with positive energy. 
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5 Conclusion 



The aim of the present study is the construction of a consistent model of a hybrid chiral bag, 
in which the exact equivalence between fermionic and bosonic (meson) phase is not assumed. 
Our results show that such a model can be actually formulated in a quite consistent manner, 
and to certain extent could be more effective way of description of low-energy hadron physics 
compared to the traditional HCM. 

First, let us note, that the initial formulation of the model is a local field theory, and 
in spite of the variety of classical solutions one needs to deal with, the covariance is broken 
only spontaneously, and so can be restored by means of the methods of refs. [23] using 
the covariant group center-of-mass variables for a localized quantum-field system. Besides 
this, the positive points of this approach are: the more correct derivation of chiral boundary 
conditions, by which any term in the initial Lagrangian has exact physical meaning; the pres- 
ence of an intermediate phase describing quasifree massive "constituent" quarks; physically 
acceptable behavior of the total energy of the bag as a function of it's geometry. Moreover, 
in this model the condition of fermion confinement, incorporated into it from the very begin- 
ning, shows up more explicitly. It manifests, in particular, in the fact that there is no need 
in the term with vacuum pressure B, which in the standard approach is inserted into the 
model by means of some extra assumptions, since in our case the Dirac's sea polarization 
itself produces the infinite increase of energy at large distances. 

The important question of the choice of method of calculation of the Dirac's sea averages 
for fermion bags should be also emphasized. The method we used is based on the discreteness 
of the fermion energy spectrum, what by means of quite obvious considerations leads to very 
simple solution of the self-consistent equations of the bag in the intermediate region. Note, 
however, that in spite of arguments in favor of such method of calculations, we cannot 
completely reject alternative methods like thermal regularization. The question of which 
one is more adequate to the physics of the problem should be answered only by means of 
detailed study of realistic models. 

Compared to ref.[21], the specific feature of the considered model is that the condition 
of nonvanishing density of the boson condensate in the interior region affects crucially the 
amount of possible bag configurations, which provide the local minima of the energy. In 
ref. [21] the boson field vanishes in the inner region, what leads to an infinite set of such 
configurations with the same topology and infinitely increasing size and energy with the main 
difference between them being the value of gradient A of the boson field in the intermediate 
region. In the present case, there exists not more than one such a configuration, and the 
excited bag states of greater energy can be obtained only by adding a number of valence 
fermions. So within the framework of 3-phase modification of the hybrid model there are 
possible quite different versions of description of composite particles like hadrons and their 
excitations. 
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7 Figures 

Fig. 1 The total energy of the bag £ft ag (a, (3) asM a function of dimensionless bag parameters 
a,P for /i = 0.25M and g = 1. 

Fig. 2 The curves of equal energy for £b a9 (a, /3)M for /z = 0.25. The minimum is clearly seen 
for a ~ 40M and f3 ~ 0.34. The values of energy in the minimum areM £b ag = 2.41, = 
1.96 and S v = 0.45. 

Fig. 3 The scheme of fermionic levels for the groundM state bag configuration. 
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